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Abstract 

We report on the Scanning Tunneling Microscope (STM) observation of linear fringes together 
with spatially varying super-lattice structures on (0001) graphite (HOPG) surface. The structure, 
present in a region of a layer bounded by two straight carbon fibers, varies from a hexagonal 
lattice of 6nm periodicity to nearly a square lattice of 13nm periodicity. It then changes into a 
one-dimensional (1-D) fringe-like pattern before relaxing into a pattern-free region. We attribute 
this surface structure to a shear strain giving rise to a spatially varying rotation of the affected 
graphite layer relative to the bulk substrate. We propose a simple method to understand these 
moire patterns by looking at the fixed and rotated lattices in the Fourier transformed k-space. 
Using this approach we can reproduce the spatially varying 2-D lattice as well as the 1-D fringes 
by simulation. The 1-D fringes are found to result from a particular spatial dependence of the 
rotation angle. 

PACS numbers: 68.37.Ef, 73. 20. At, 73.43. Jn, 61.72.Nn 
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INTRODUCTION 



Since the discovery of the scanning tunnehng microscope (STM) by Binnig and Rohrer 
highly oriented pyrolytic graphite (HOPG) has been one of the most commonly used 
model surfaces for investigating and modeling of the STM imaging process. HOPG has a 
hexagonal structure with weakly coupled layers of graphene stacked in ABAB... sequence. 
Graphene is a semi-metal with zero density of states (DOS) at Ep but a weak inter-layer 
interaction in graphite makes this DOS finite but small. This small DOS is quite sensitive 
to defects that influence the overlap of orbitals between layers [2|]. This makes such defects 
easily visible by STM. Further, the weak inter-layer interaction makes such defects, like 
moire patterns, possible near the surface of HOPG. A widespread study of this material has 
led to the discovery of several types of crystal imperfections native to the basal plane of 
30PG. These include stacking faults, graphite strands and fibers, broken or flaked layers 
I .upe.la«ces MM with pe„od.city of about ten _te. o. even la.ge, and 
the buckling of the top layer [8|]. 

Broadly speaking, super-lattice structures in HOPG can occur due to three reasons. First 
type appears in graphite-intercalated compounds [9| with fixed period super-lattices like. 



or 2x2. Second type is the electron density oscillations (h 



Id, 



ie Friedal oscillations) 



These oscillations 



near point defects with a periodicity times that of graphite 
decay rapidly as one moves away from the defect. The third one involves disorientation 
of one or more of the graphite layers near the surface [sl, y, [?]. This type of super-lattice 
structure is commonly known as moire pattern and it arises as an interference pattern of two 
identical but slightly rotated periodic lattices. As first suggested by Kuwabara et. al. j^, 
this rotation by a small angle 9 gives rise to a super-lattice of the same symmetry but with 
much larger lattice constant given by, D = d/(2sm{6/2)), with d as the lattice constant of 
the original lattice. This super-lattice is rotated with respect to the original lattice by an 
angle = 30° - 9/2. 

In this paper we present STM study of a spatially varying 2-D super-lattice pattern 
and 1-D fringes observed on HOPG (0001) surface near some defects. Such 2-D patterns 
have been reported earlier by several groups and are widely believed to be due to the moire 
interference 

mm 

. Spatially varying periodicity patterns have also been reported before 
and have been attributed to a shear strain [12] or quantum confinement in a linear potential 
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[l3|. The pattern observed here seems to be a result of an in-plane shear strain causing a 
spatially varying moire rotation of a top layer. However, this is the first observation of linear 
fringes connected with a moire pattern. We also discuss how these fringes can arise from 
a shear strain using a simple theory similar to the one proposed by Amidror and Hersch 
[isl for moire patterns in optics. This approach gives us a better analytical insight into 
the large scale structure of such patterns. An earlier model on moire patterns in graphite is 
purely computational, based on the variation in local density of atoms (lo}]. In our model, 
the moire patterns are examined in Fourier transformed k-space to quantify the variation 
in local stacking for a fixed as well as a spatially varying moire rotation angle. Using this 
idea we can understand and simulate the above observed pattern and, in particular, the 1-D 
fringes. 

EXPERIMENTAL DETAILS 



Experiments were done with a home built STM similar to the one described elsewhere 

I]. The data reported here 



17l |. This STM uses a commercial electronics and software [18}] 
were taken in ambient conditions. HOPG was fixed on the sample holder with a conducting 
epoxy and the sample was freshly cleaved using an adhesive tape before mounting it on 
the STM. Fresh cut Pto.8lro.2 wire of 0.25 mm diameter was used as the STM tip. The 
images have been obtained in constant current (feedback on) mode. The images shown here 
are filtered to remove steps and spikes, however, for quantitative analysis we have used the 
unfiltered data. 



RESULTS 

Spatially varying lattice and 1-D fringes 

FigilK shows a topographic image of (0001) basal plane of HOPG with the super-lattice 
structure together with two steps and two carbon fibers. We mark these as fiber-1 & 2 and 
step-1 & 2. The step-1 is a double layer step with a height of 0.7±0.1nm while the step-2 
is a single layer step. These step heights have been found from the images taken away from 
the super-lattice pattern area. These two steps divide the area into three terraces with the 
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FIG. 1: a. A STM image (0.32x0.32 fim'^ at 0.5V and 0.2nA) showing a spatially varying super- 
lattice structure confined between fiber- 1 and fiber-2 on the layer defined by step-1 and step-2. 
Near the termination of superstructure a 1-D fringe pattern appears, b. A later STM image 
(2.16x2.16 firri^ at 0.8V and 0.2nA) showing a larger portion of fiber-1 and fiber-2. The area of 
left image marked by the white square. 

middle one having a super-lattice pattern confined between the two fibers as seen in the 
same figure. 

Carbon fibers have been observed earlier by various groups on the graphite surface jj]. 
These, we believe, are long and thin ribbons or rolls of graphite created, presumably, during 
the cleaving process. The exact vertical location of the fibers relative to various visible 
layers is not so clear in this image. However, since the fibers are strictly limiting the pattern 
on the mono-layer terrace defined by step-1 and step-2, we believe that these fibers are in 
contact with this affected layer. They may be either below or above this particular layer. 
If a layer is above the fiber it would go over the fiber with some buckling. Such buckling 
cannot be ruled out by stress considerations as the fibers are quite wide (>10nm) and not 
so high (<lnm). 

A relatively larger area image including the patterned region is shown in Fig{T]D. This 
image was taken later when the fiber-1 was fully relaxed as discussed in detail later. The 
region of the left image is marked by a square in the right one. The other ends of the two 
fibers are also visible in FigHb. Here we see that the fiber-2 looks brighter above the step-1 
and the fiber-1 has similar behavior but this is seen in the FigH^. Here, both the fibers 
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FIG. 2: The topographic images of the super-lattice structure showing its evolution with the motion 
of fiber-1. a. 0. 52x0. 52 fim^ , b. 0.35x0.35/xm2, c. 0.43x0.43^m2 taken at 0.5V and 0.2nA. 

look less bright when they emerge from the covered layer(s) near the center of the image in 
FigHb. 

As shown in Fig{T^ the super-lattice contains two regions namely a 2-D lattice extending 
till the end of fiber-1 and a 1-D wave-like pattern starting from the fiber-1 end. The 1- 
D fringes bend towards the step-2 and terminate at this step. The 2-D lattice continues 
beyond the fiber-2 end but terminates at step-2 in this image. The 1-D pattern contains 
the same number of maxima and minima as the terminating 2-D lattice. The 2-D lattice is 
not uniform as it evolves from a hexagonal lattice deep inside the fibers to nearly a square 
lattice with much larger periodicity. The super-lattice structure is also present over the 
step-1 (up to fiber-1) but with corrugation reduced by a factor of 2.3. Although this is not 
clearly visible in FiglUbut we clearly see it on zooming into this area. The corrugation is 
also found to vary with periodicity as analyzed later. 

Pattern variation with fiber-1 motion 

As we image this area repeatedly, fiber-1 becomes shorter in length as the whole fiber 
seems to recede. This may be due to some stress on the fiber. Several groups have reported 
various surface modifications in the STM experiments but the exact mechanism behind 
such modifications is not quite understood. We show three of the latter images in Figl2] 
depicting the changing pattern with the fiber-1 withdrawal. The withdrawal of this fiber 
is not uniform with time as sometimes a significant length disappears in one scan and at 
other times it recedes gradually. Eventually the receding stopped and the length of the 



fiber stayed same over a couple of days. One remarkable observation is the pinning of the 
boundary between the 2-D super- lattice and the 1-D fringes to the fiber- 1 end. With the 
receding of fiber-1, the 1-D fringes terminate on the fiber-2 as opposed to step-2 in earlier 
images. Moreover, the spatial extent and corrugation amplitude of 1-D fringes increases 
together with some modifications in the large scale structure of the 1-D pattern. 

ANALYSIS 

Two close-up images of this pattern are shown in FigI3] in two different regions of the 
structure at different orientation of the image window as compared to the previous ones. 
These images have been taken at an early stage when fiber-1 was extending the most. As 
apparent from the image shown in Figj3^, the pattern deep inside the fibers (towards left) is 
a hexagonal lattice of fixed periodicity. As we move out, the periodicity increases with the 
pattern evolving to an oblique lattice and then to nearly a square lattice as seen in Figj3b. 
Fourier transforms (FT) of these images as shown in the inset also illustrate this variation. 
In detail, the FT in FigJS^ shows six sharp spots while the FT in FigjSjD has these spots 
elongated. 

To analyze the periodicity variation carefully, line cuts were taken along various lines 
perpendicular to the line 'C in FiglSb- Two such representative lines are A and B along 
which the topographic height is plotted in FigOb and FigOli. These lines were carefully 
chosen to pass over the bright spots of the 2-D lattice. The periodicity along the lines 
perpendicular to C was found to be visibly constant but as one moves out along C this 
periodicity increases. The corrugation along the line C also varies. The RMS corrugation 
and periodicity as measured by line cuts perpendicular to C are plotted in FiglS^ and FigjSf 
as a function of distance along the line C. The periodicity along such lines has been found 
from the average separation between peaks. 

If we use the moire rotation hypothesis, the rotation angle 6 (between the top layer and 
the layers underneath) should also change spatially to get a varying lattice spacing D. Using 
D = d/{2 sin(^^/2)) (d=0.246nm), we plot the variation of 6 in Fig|3^. This varying 9 implies 
that this rotated layer undergoes a shear as we move out of the fibers. The portion of this 
layer that is well inside the fibers has a constant rotation and the portion far outside the 
fibers is free of any moire pattern and thus 6 is zero. The intermediate region is strained 
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FIG. 3: a. Super-lattice pattern image (120xl20nm^ at 0.5V and 0.2nA) taken deep inside the 
fibers, and b. image (196xl96nm2, 0.5V/0.2nA) taken near the boundary of the pattern. The 
insets show the respective Fourier transforms, c and d show the line scans along the lines A 
and B. e, f and g show the variation (along line C) in rms corrugation, periodicity, and rotation, 
respectively. 

with a changing 9 and periodicity. As estimated from the periodicity, the angle 9 changes 
from 2.3° to 1.3°. The 1-D fringe pattern which has a close correlation with the 2-D moire 
pattern prompts us towards a common origin for the two patterns. We believe that ^'s 
particular spatial dependence is responsible for the 1-D fringes as discussed later. 

A second order polynomial fit (as shown in Figl3^) for this 9 variation gives an expression, 

9 = Aq-Aix- A2X^ (1) 

with = (2.11 ± 0.02)degree, = (3.0 ± 0.6) x lO'^degree/nm, ^ = (2.4 ± 0.5) x 

1 "^degree /nm^. This is also the expression for local angular orientation for a cantilever 

□ 

with one end clamped and the other end loaded [19|. 

A possibility on how the stress arises in the graphite layer is similar to the bending of a 
cantilever with certain loading conditions caused by the fibers. By the definition of strain, 
the maximum strain at a given x would be To estimate the maximum strain, the ^ 
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is maximum at the boundary of the 2-D super-lattice and 1-D fringes and it is found from 
the polynomial fit (EqU]) to be (1.6±0.3)xl0^^rad/nm. At this location the thickness (d) 
is 115nm. This point gives the largest stress, crmax=l-l±0.2 GPa, using F=121.9 GPa [20| 
for graphite. For comparison, the macroscopic indentation experiments indicate that stress 
on the order of ~1 GPa damages the surface of graphite 2l|]. On the other hand, Snyder 
et. al. [20I have shown that the stress required to induce dislocation motion on basal plane 
of HOPG is of order 5-200 MPa. 



UNDERSTANDING THE LINEAR FRINGES FROM SPATIALLY VARYING 
MOIRE ROTATION 

Here we describe a simple model to understand and to simulate the large-scale structure 
of the spatially varying moire patterns. For illustration, let us first consider a 1-D periodic 
lattice. If we superimpose two 1-D periodic patterns given by cos kix and cos (| ki—k2 |^ 
ki, ^2)5 we would see an interference (or beats-like) pattern with periodicity 27c/\ki — k2 \ in the 
superposed pattern (i.e., cos fcix cos /c2a;). In 2-D if we superimpose two hexagonal patterns 
of slightly different reciprocal lattice vectors ki and k- (i = 1-6), the resulting pattern in 
real space would come from the difference of the two sets of lattice vectors, i.e., (k- — kj). 
The longest period variations would arise from Aki = k- — ki (see FigJlj), which are actually 
responsible for the moire patterns here. In 2-D, the pattern could arise due to differences in 
both the magnitude and direction of k and k' |lj] but in the present case only the direction 
mismatch is playing the major role. Aki as seen from FigJH has a magnitude of 2A;sin(^/2) 
and the smallest angle it makes with one of the ki's is 30°— 9/2. This is consistent with the 
aforementioned periodicity and angle of the moire super-lattice. Now the moire interference 
pattern in real space is easily produced by inverse Fourier transformation (IFT) of these six 
Aki vectors. This IFT gives a pattern described by, 

I{x,y) = Jo + /i[cos(/i) + cos(/2) + cos(/3)], (2) 

where /q and Ii are two constants and 

/„ = 2A;sin(^/2){-|/sin((n - l)7r/3 + ^/2)+xcos((n - l)n/3+9/2)}]. (3) 

Here n=l, 2 or 3 and k = Air/^/Sd with d (=0.246nm) as the real space ab-plane lattice 
parameter of graphite. Each cosine term in Eql2] gives rise to a 1-D periodic pattern with 





FIG. 4: a The lattice vectors of the fixed (k) and shghtly rotated (k') hexagonal lattices and their 
differences (Ak). In this figure y-axis is pointing vertically up and x is pointing to the right, b 
Moire interference pattern calculated using Eql2]with a 6=2.5°. The size of the image is 25x25nm^. 



bright and dark fringes along parallel straight lines (for constant 6). The bright fringes of 
cos fn are described by the contours fn{x,y) = 2iV7r, with N as an integer. For a constant 
6, the fringes due to the three cosine terms make an angle 7r/3 with each other and the 
sum of these three cosines leads to a 2-D triangular lattice pattern. One such pattern is 
shown in FigJDD for 9 = 2.5°, /o=0, and /i=l. In this model, the I{x, y) only quantifies the 
local stacking pattern. The I{x,y) has a variation between +3 and -3 (taking Iq = and 
Ii = 1) with a value of +3 signifying an AA stacking and -3 signifying an AB stacking. An 
intermediate value represents a slip with the degree of slip quantified by the magnitude of 
I{x,y). The DOS at the surface is a result of top three layers and so the DOS at BAB and 
BAG could be actually different while our model does not differentiate between these two 
stacking sequences [s,!^]. 

This model can :.s„ be used f„. tKe .eo.et.eal.y .a„sfo..ed .at.ee. fl, wKe.e the 
rotation angle is spatially varying. This way we can model the spatially varying moire 
patterns and understand their large scale structure. In particular, we can find an analytical 
condition on the spatial variation of the rotation angle that gives rise to the 1-D fringes as 
following. For spatially varying 6, the 1-D fringes due to each cos /„ will have a spatially 
varying periodicity and, in general, these fringes may not be straight and parallel to each 
other. These fringes will be locally perpendicular to V/n and their local periodicity can be 
quantified by 27r/|V/rt|. These /„'s (see Eqj3]) satisfy the condition, 



/i — /2 — /s- 



(4) 
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For I{x, y) to have a 1-D pattern, the fringes due to the contributing cosine terms must 
be locally parallel to each other or, in other words, each /„ should have the same constant 
value contours. In particular, and as it turns out to be the case for the observed 1-D fringes, 
if the 6 variation is such that /i = C (a constant) then using EqJU /2 = C + fs- For 
such 6 variations cos /i gives no fringes while the fringes due to cos /2 and cos /a are locally 
parallel with the same local wavelength. Thus I{x, y) would have only 1-D fringes. Similar 
argument would also hold for /2=constant or /3=constant, which is obvious from symmetry 
considerations. The expression for fi (Eql3]) can be simphfied to. 



In small 6 approximation we can neglect the {cos 6 — 1) term in comparison to sin 6' term in 
Eq|5]and this gives 6{x,y) = C/kx. Here, one has to be careful in choosing the origin while 
using this form of 6{x, y). The same origin has to be used in calculating the J(x, y) to ensure 
1-D fringes. In fact this choice of origin gives one more parameter in choosing 6{x,y). 

A simulated pattern using above ideas is shown in FigJ^b together with the observed 
pattern in Fig|5^. Here, we have used two different spatial forms of 6 in different x intervals. 
We have used 6{x,y) as found in Eq{T] for 2-D region, where 6 changes from 2° to 1.29°. 
We use the small 6 approximation for 1-D fringes, i.e., d{x,y) = C/k{x + a) for 6 between 
1.29° and 0.76°. The constants C/k and a are found to be 3.067 rad-nm and 4.399 nm, 
respectively, by ensuring the continuity of 6 and dd/dx at x=132.3 nm. This is necessary to 
keep the stress finite and continuous everywhere. For 1-D fringes the required 6{x,y) is of 
the form of C/kx, so in FigJSl we have actually shown images of I{x + a, y) with a = 4.399 
nm. A plot of the used 6 variation as a function of x is shown in FiglSt. Since the 2-D to 
1-D pattern boundary is pinned to the fiber-1 end we believe that there is a sudden change 
in spatial dependence of stress across this x causing a change in 6 behavior that in turn is 
responsible for the boundary. FigO also shows the three cos /„ maps, where one can see 
clearly how the pattern arises in the two regions in terms of the three cosines as argued 
before. In the 1-D fringe region, we see that cos /i gives no fringes while cos /2 and cos fs 
have fringes that are locally parallel to each other. In the 2-D pattern region the fringes due 
to the three cosines are not parallel. 



/i = kxsinO + ky{cos6 — 1) 



(5) 



and so for y) = C v/e find. 
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FIG. 5: a. STM image of the super-lattice pattern of varying periodicity and linear fringes 
(196xl39nm^ at 0.5V and 0.2nA). The 2-D pattern is present up to 100 nm and beyond this it 
changes into linear fringes b. Simulated image (see text) of size 198xl41nm^ with the 6 variation 
shown in c. d, e, f show the three cosine components, i.e., cos/i, cos/2 and cos/3, respectively, of 
the image in b. (In figures b, d, e and f y-axis is pointing at angle 3° to vertically up and x-axis 
is perpendicular to y and towards the right. 

So far we discussed a particular condition on 6 that gives 1-D fringes. Now we find a more 
general condition on 6 for getting 1-D fringes and we also argue here that above used spatial 
dependence of 6 for 1-D fringes is the only possibility (in small 6 approximation) for the 
observed fringes. The existence of 1-D fringes requires that all the constant value contours 
of /i, /2, and /s are identical. This is possible if fi = ^(/a) and, using EqJU /2 = /s + 0(/3), 
where is a scalar function. If we assume to be a linear function, i.e., (f){x) = a{x + Ci) 
then /i — af3=aCi. This, using Eq|3l gives the same solution for 0{x,y) as Eq|6] except 
for X, y, and C replaced by x' = x{a-^ + 1/2) - y{V5/2), y' = y{a~^ + 1/2) + x(V3/2) 
and Ci, respectively. In small 6 approximation, this gives 6{x,y) = Ci/kx'. Therefore, any 
function of the form 6{x,y) oc l/(x + 7?/) would give rise to linear fringes for arbitrary 7 
(=\^a/ {a + 2)). The form of 9 used in producing the linear fringes in Figl5] is a special 
case with a=0 and Ci ^ 00 with aCi=C. From Fig|5K we see that the observed 1-D fringes 
connect smoothly to the rows of bright spots coming from cos /2 and cos /s while the rows 
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of bright spots due to cos /i terminate before the 1-D fringes. This means that cos /i is not 
contributing to the observed 1-D fringes and cos/i has to be constant in this region. So 
6 = C/kx is the only possibihty for the observed fringes in small 6 approximation. 

DISCUSSION 

Our simulation nicely captures most of the details of the observed pattern except for the 
curvature of the 1-D fringes. For the curvature in 1-D fringes, we find that any modification 
in spatial dependence of 6 from C/kx affects the 1-D nature more seriously than producing 
the curvature. Incidently, using the exact form of 9 as in Eq|6]as opposed to its approximate 
form of C/kx does not give the curvature either. Here, we believe that our approach of 
varying only the directions of k's and not their magnitude is a bit oversimplified though it 
captures the general structure of both the 1-D and 2-D patterns. However, to get both the 
direction and magnitude of k's requires a detailed understanding of the strain field in the 
affected region. This is not possible from the STM images alone as we have to know the 
stresses on the boundaries to solve the complete boundary value problem using elasticity 
theory. 

As pointed out earlier, we should be cautious in using this model as this does not describe 
the quantitative contrast of an STM image. For instance, this model cannot explain the 
variation in corrugation amplitude of the 2-D super-lattice with its periodicity as seen from 
Figl3^. The STM contrast here actually represents a DOS variation as our local tunneling 
spectra show that the bright regions are more metallic than the dark ones consistent with 
earlier work by Kuwabara et. al. [5]. As of now no calculations on electronic DOS on 
moi re p atterns in graphite exist. A computational model was proposed by Hentschke et. 
al. [23] and reviewed by Pong et. al. [16], which is based on the variation in local density 
of atoms. In a sense our I{x,y) also quantifies the local density of atoms but our model is 
more appropriate for the analytical understanding of the large scale structure of the spatially 
varying moire patterns. 
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CONCLUSIONS 



In conclusion, we have studied a spatially varying 2-D and a connected 1-D super-lattice 
structure on HOPG. This is the first observation of such 1-D fringes connected to a 2-D 
moire lattice and we attribute this pattern to a spatially varying moire rotation of a top 
graphite layer. The spatially varying rotation implies a shear strain in this layer. We have 
also described a simple model for understanding spatially varying moiree patterns. This 
model can successfully simulate the observed pattern by using a spatially varying rotation 
angle. The 1-D fringes are found to arise from a particular spatial dependence of the moiree 
rotation angle. 
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